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Abstract 

I compute non-perturbative corrections to the kernel governing the 
evolution of non-singlet parton densities. The model used is QED in 
the limit of many charged particles. I find an ultraviolet renormalon 
corresponding to a term of order Q 2 /A 2 , where A is the pole of the 
coupling constant. This term has a non-trivial dependence on the 
variable x = —q 2 /(2pq) and its coefficient scales as rr 3 /(l — x) 2 (p is 
the momentum of the hadron and q is the momentum transfer). An 
extrapolation of my results to QCD implies a breakdown of the parton 
model near the elastic region. 
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1 Introduction 



One of the main processes for the study of the hadron structure is deep 
inelastic scattering: 

I + H -> I + X (1) 

where if is a hadron, I is a lepton and X is any hadronic final state. 
In the parton model picture, the hadron is viewed as a gas of fast moving, 
real, collinear and non-interacting partons (quarks and gluons), which share 
the energy-momentum of the hadron. Perturbative QCD accounts for hard 
gluon radiation and produces logarithmic corrections to the naive parton 
model. The cross section a can be written as the product of a short-distance 
('hard') cross section times a parton density: 

a(p, q) = a hard (xp, q) ■ q(x, Q 2 ) (2) 

where Q 2 = —q 2 > 0. 

In perturbative QCD one can compute the hard cross section and the evo- 
lution of the parton densities from a scale Q 2 to a different scale Q' 2 |], @]. 
The parton densities at a given scale cannot be predicted and have to be 
determined experimentally. 

There are corrections to the parton model which originate from the in- 
trinsic transverse momenta of the partons, their virtualities, 

as well as from parton-parton interactions. These phenomena are all 
characterized by the hadronic scale Aqcd ~ 300 MeV. A kinematical anal- 
ysis shows that taking these effects into account produces power-suppressed 
corrections to the factorized cross section (|2]), of the form H 

where n = 1, 2, 3, k, ... (3) 

It is crucial to have an estimate of these non-perturbative terms for a trust- 
worth application of perturbative QCD. 

The occurence of power suppressed corrections can also be demonstrated 
with the operator product expansion (OPE) 0, [|. The cross-section of the 
process ([]]) is proportional to the imaginary part of the forward scattering 
amplitude 

(H | V I H) (4) 
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where 

h^ u = ij tfxe-*** T Ux)Jl(0) (5) 

is a hadronic tensor and J M is an electromagnetic or a weak current. 
The tensor h^ v can be expanded in the limit, 

Q 2 , p-q — ► oo, a; = constant, (6) 

in a tower of covariants of progressively higher twist T v(T = dimension of 
the operator - spin): 

Ku = £a(Q7^)c^ 2 ) (7) 

where Cj(Q 2 //i 2 ) are coefficient functions, Oj(/z 2 ) are local operators and // 

is a subtraction point. 

The lowest twist operators (T = 2) 

i/)'yp\aD 1 ,D lil D tt2 ...D lin ip - traces, F^D^D^.-.D^F^ - traces (8) 

correspond to the parton model, while higher twist operators produce cor- 
rections of order 1/Q 2 (A a is a matrix in flavor space). 

Arguments based on parton kinematics or the OPE can provide only order 
of magnitude estimates of the terms @. It is the purpose of this paper to give 
a semiquantitative estimate of the power suppressed corrections, based on the 
renormalon calculus. In general, there are power suppressed corrections to 
the evolution kernel and to the hard cross section and the observable effect 
is a convolution of both terms. We work hovewer in the leading logarith- 
mic approximation (see eq. p0|)), in which the power suppressed corrections 
to the hard cross section are absent. The power suppressed corrections to 
the evolution kernel are instead non-vanishing and have a direct physical 
meaning. 

The idea of the renormalon approach is that the complete perturbative 
series reflects non-perturbative properties of the process |J. The chain of 
arguments is the following. The perturbative series is strongly divergent in 
all the known cases, i.e. it diverges for every non-zero value of the coupling 
constant. The coefficients have typically a factorial growth, 

\cn | ~ nl (9) 
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A meaning can be assigned to the perturbative series through a resummation 
of the series itself. This operation requires in general an arbitrary prescrip- 
tion, which introduces an ambiguity. The latter is a size of the intrinsic 
non-perturbative corrections. 

Let us present a physical argument of why the investigation of the power sup- 
pressed corrections (^j requires the knowledge of the complete perturbative 
series. Consider a quantity that has a simple expansion in powers of a s (Q 2 ) 
(i.e. an infrared safe quantity). Since 



Therefore, it is consistent to look at the power suppressed corrections only 
after the inclusion of all the perturbative terms. The renormalon calculus 
cannot be applied to QCD because we are not able to compute the coefficients 
of any observable for any n. The use of models cannot be avoided and various 
kinds of 1/N expansions have been used. We consider a 1/N expansion 
of QED, where iV is the number of charged particles 0. This model is 
indeed solvable and involves tree diagrams containing strings of self-enegy 
corrections to the photon lines (see fig. 1). Since QED contains massless 
vector quanta like QCD, collinear singularities are the same in both theories. 
Therefore, perturbative physics is the same as in QCD as far as charge flow 
effects are neglected. 

Actually, QCD can be modelled by changing the sign of the beta function 
at the end of our QED-like computation || . In such a way the infrared prop- 
erties of QCD are characterized by the perturbative growth of the coupling 
constant at small momentum transfer. Loosely speaking, we are dealing with 
tree diagrams in which the bare coupling is replaced by the QCD one-loop 
running coupling constant. 

The paper is organized as follows. Section (§) contains a short review of 
the Borel transform technique and of the renormalon calculus. In section (^) 
we compute the input diagrams for the kernel P qq . The latter is derived in 
section ([4]). Section (|5|) contains a discussion of the results and section (^) 




(10) 




for any positive n and m 



(11) 
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traces the conclusions. There are also two appendices. In appendix A the 
technique to compute the integrals of the Borel transform is described, while 
in appendix B we prove some properties of the planar gauge. 

2 The Borel transform 

The perturbative expansion of an observable O in quantum field theory is 
given by a series of the form 

oo 

0(a) = c+Y,c n a n+1 . (12) 

n=0 

This series is strongly divergent in all known cases because the coefficients 
have a factorial growth with n: \ c n |~ n\. 

As a O-dimensional model of quantum field theory, we can consider the inte- 
gral: 

I k (a) = / dx x k e~ x ~ ax (13) 

J —oo 

There is an instability for a < and the expansion in powers of a is divergent. 

By explicit computation one finds c n = (— l) n+1 T(2n+fc/2+5/2)/r(n+2). 
For large n there is a factorial-like growth of the coefficients with alternating 

sign: c n ~ (_l)»+l2 2 "+ fc /2+3/2 e (fc-l)/2( n + !)|. 

The simplest interpretation of the expansion (O) is that of an asymptotic 
series. It is sensitive to truncate the series at its smallest term, which is 
a measure of the error. A more detailed analysis aims to reconstruct the 
original function from the knowledge of the c n 's. This is achieved by means 
of the Borel transform technique. The Borel transform is defined by: 

6{t) = c 8(t) + t°^t n (14) 

71=0 U - 

We assume that the series on the right-hand side defines an analytic function 
of t in a neighborhood of t = 0, which can be analitically continued in all the 
Borel plane. The observable O can be reconstructed from the Borel trasform 
by means of the formula: 

B (a) = / dt e~ t/a 0(t) (15) 
Jo 
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If, for example, c n = (— l) n n\, the series ( |T2| ) has zero radius of convergence, 
while the Borel transform defines an analytic function in all the t-plane, given 
by 0(t) = 1/ (1 + t). Since 0(t) has no singularities for t > 0, the observable 
can be reconstructed unambigously by means of the inversion formula (|15|). 
In such a case we say that the perturbative series is Borel summable. On the 
contrary, if c n = n\, the Borel transform is given by O (t) = 1/(1 —t) and there 
is a singularity at t — 1. The integral in eq. (|H|) must be supplemented with 
a prescription to deal with the pole. The arbitrareness of the prescription 
reflects the intrinsic ambiguity in the reconstruction of 0(a). 

Let us discuss now how renormalons are related to power suppressed 
corrections in quantum field theory. Consider a Borel transform of the form: 

K 

0(t) = + regular function. (16) 

The ambiguity in the observable is given by: 

50 = K e~ nMQ2) (17) 

where we have taken the pole residue as a measure of the ambiguity. 
Expressing the coupling constant in terms of the relevant scales (dimensional 
transmutation), 

a{Q2) = -Mcgy/A" (18) 

we have: 

/ A 2\ -bon 

50 = K — , (19) 



where bo is the first coefficient of the /^-function, defined including the factor 
1/(4tt). In QED, b = 1/(3tt). 

The position of the renormalon pole determines the power of the non- 
perturbative correction, while the pole residue sets the coefficient of the cor- 
rection. In general, the Borel transform has a collection of singularities along 
the real axis and the leading ambiguity is related to the pole on the positive 
axis closest to the origin. 

Up to now, the existence of renormalons has not been proved in any realistic 
theory (such as QED, QCD, A0 4 , etc.). We consider QED at order 1/N. 
In this model the existence of renormalons is related to the Landau pole, as 
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shown by the following intuitive argument. The amplitudes are given by loop 
integrals containing the one-loop running coupling constant (18) 



"•°<*>g (* + ^M* + V " (20) 

where Pi,P2 ■ ■ ■ denote external momenta. The coupling constant diverges at 
the Landau pole, producing a divergence of the integral. The latter has to 
be regularized, and this operation introduces ambiguities in the amplitude. 
There is also a more explicit way of understanding the relation between 
Landau pole and renormalons. Expanding the running coupling constant in 
powers of a at a fixed scale /i, 

a(k>) = - ^ 



1 — boa(n 2 ) \og(k 2 /n 2 ) 

oo i 2 

E (b \og-Ta(v 2 ) n+ \ (21) 

n=0 Z 1 



gives rise to a coefficient of a(/i ) n of the form 



1 dx flog-) = n\ (22) 



V X 



There is a fixed sign factorial divergence resulting in renormalons in the Borel 
plane. 



3 The computation 

In this section we evaluate the input diagrams for the derivation of the evo- 
lution kernel P qq . The model is QED with N flavors, and we work at order 
1/N. Quarks are assumed massless and the incoming quark is taken on-shell, 
i.e. p 2 = 0. Soft singularities (infrared and collinear) are regulated by the 
Borel variable t itself. For t > the gluon propagator is indeed less singular 
than in the usual case (see appendix A for a formal proof of the regulator). 
It is convenient to use the planar gauge with the gauge vector directed along 
the final quark direction. In this gauge collinear singularities decouple from 
the final quark leg and the only relevant diagram are the ladder ones (see 
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appendix B for a proof). The free photon propagator is given by: 



'vm - K (23) 

where 

S^k) = -9„+ Knu+ ^ n " , (24) 

is the polarization sum, and n = xp + q is a light-like vector, n 2 = 0. 

The lowest order diagram is given in fig. 2, and is of order (1/N)°. The 
rate is given by: 

W ° = 2ff(p)Q 2 ^ 1 ~ ^fcPyvixp + g)7„] (25) 
where -E(p) is the energy of the incoming quark. 

The next order diagrams involve a single photon exchange, and are there- 
fore of order a = a/N. Inserting bubbles into the photon line involves the 
coupling a = aN. Diagrams with multiple gluon emissions are of order 
a k = a k /N k with k > 2. We work at order 1/N and to all orders in a. 
This implies that we have to consider the diagram of fig. 3, with the photon 
propagators containing strings of an arbitrary number of bubbles. 
The proper self-energy at one-loop is given by: 

m^(k) = i(k^k v - gflu k 2 )n(k 2 ) (26) 

where 

U(k 2 ) = -ab log(^p^ (27) 

in the momentum subtraction scheme. 

The full photon propagator is therefore given by: 

iA^fc) = + mo°(k) ^® + ... 

k 1 + ie k 2 + le k z + le 



[k 2 + ie] [1 + Tl(k 2 )} 



(28) 



The inclusion of the bubbles modifies the free propagator by the factor 1/[1 + 

n((p-*) 2 )]. 
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The diagrams involved at order 1/N are given in flg.4. For diagram (a) 
we have: 



M. = e u(k + q y,^ %u{p) i ^-_ k l )2) ^ - k) (29) 



The contribution of this diagram to the transition probability per unit time 
is given by: 

Wa = 4nNE(p) I ^ k Tr [H>P7<xHXfc + Vhn] x 

x f(a) 5 + {{k + qf) 6+{(p - k) 2 )S pa (p - k) (30) 

where 5 + (x 2 ) = 5(x 2 )9(xo) and we defined 

/(a) = |i + n((p-^) r (31) 

The dependence on a is contained in this function, which can be expanded 



as 



/M - | |^W(-^)"(Mo g ^) a—. (32) 



According to the definition in eq. fll4|) , the Borel trasform is given by: 

/(i) .( cos .„_^ log (^)(_Z_)" ,33, 

where u = —b t (notice a difference of sign with respect to the usual defini- 
tion). 

The Borel transform w a is computed substituting f(t) in eq. (|30| ) in the place 
of f(a). Because of the S + ((p — k) 2 ), w a vanishes for u < 0. By analytic 
continuation (we assume that the Borel transform is an analytic function), 
the Borel tranform vanishes in all the Borel plane: 

w a = (34) 

This result is somewhat paradoxical because it implies the vanishing of the 
contribution of the basic ladder diagram, which is known to produce the one- 
loop scaling violations in QCD. There exists however a physical explanation 
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to eq. (0). The dynamics in our model is controlled by the distribution of 
invariant masses (p — k) 2 flowing into the fermion bubbles. In diagram (a), 
we have that (p — k) 2 = 0, implying that any non-trivial contribution to the 
Borel transform is impossible. 
The amplitude of diagram (b) is given by: 

2 / t i \ k r \ Spcr(p k) 



k 2 + ie IP ^ J [(p - k) 2 + ze][l + U((p - k) 2 )] 

u(p-k- l)jaV(l) (35) 



The partial rate is given by: 

1 f , 4 , S»°{p-k) 



12tc 2 NE(p) J (p- kf 

x TrP + q) ltl k lp p lc k lv ] 5 + {{k + g) 2 ) 0+((p - A;) 2 ) (36) 

where 6 + {x 2 ) = 9{x 2 )9{xq) and we defined 

9{a) = \ i+mp-k) 2 ) \ 2 (37) 

The phase-space integral over the final quark-antiquark pair has already been 
done. 

The dependence of w b on the coupling constant a is in 

»«0 = | !^f (-W(Mog^)V- (38) 
whose Borel transform is given by 0: 

To compute w b , we substitute the expression of g(t) in eq. (|36|) in the place 
of ff(a). 

In the framework of a 1/iV expansion, the leading log approximation 
(LLA) involves the resummation of all the terms of the form: 



l , Q 

Trlog- 
N p 



log^ , fc = 0,l,2,3..,n,... (40) 
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In the evaluation of the kernel, only those terms in w b which diverge loga- 
rithmically at u = have to be considered. This means that we need only 
to keep the terms in the trace proportional to k 2 or to k ■ p. It is easy to 
see that k ■ p is of the same order as k 2 . The collinear singularity is indeed 
produced by the virtual states of the photon with a small invariant mass, 
i.e. (p — k) 2 ~ 0, which implies 2p ■ k ~ k 2 (see appendix A for a rigorous 
proof). The terms proportional to (k 2 ) 2 , k 2 k ■ p, etc. do not give rise to 
any logarithmic singularity and therefore do not contribute to the evolution 
kernel (they contribute to the 0(1/N) hard cross section). Therefore, once 
a power of k 2 or k ■ p has been extracted, the trace can be evaluated in the 
collinear limit, k 2 = 0. This condition together with {k + q) 2 = 0, implies 



A: 



xp. 



The rate is given by: 



2u 



W b 



sin 7TU /i 

67r 3 b NE(p) 



1 O 2 

X X 



Tr[xp~f u (xp + g)7 P 



(41) 



(42) 



where Q 2 = —q 2 > 0, and the integrals ii, I 2 and J 3 are given by: 



h 



-I 
-I 
-I 



d 4 A; 



d 4 A; 



d 4 A; 



5 + {{k + qf) 6 + ((p-kf 



6 + ((k + q y 



2\2 



(k 

5 + ((k + qf 
k 2 



2p ■ k 



(( P -k) 2 y+ u 

6 + ({p-k) 2 ) 
((p - k) 2 ) l+u 

6 + {{p-k) 2 ) 
{{p-k) 2 ) 1+u n-(p-k) 



(43) 



The final result is (the computation of the integrals is described in detail in 
appendix A): 



w b = 
where 
A(x, u) 



sin iru fj, 



2u 



l2ir 2 b NE(p) 



Tr[xpy v (xp + q)n/»] 



x 



u+l 



(Q 2 Y+\l-x) 



— A(x,u) (44) 



lit' 



+ 

xu 



+ 2 



logx 



1) 



1 + u - - J F(l, 1 - u, 2 - u; 1 - x) + 



■F(l,-«,1 



+ 



1 — x 



1 — x) 2 u 
- LerchPhi(— x, 2, — u) + 



u: 



-x) + 
2 



l-x) 



h(x, u). 



(45) 
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The special functions are defined in appendix A. 

Let us end this section with some comment about the decay rate (44). The 
single pole at u = is the collinear singularity and replaces, loosely speaking, 
the pole in e of usual dimensional regularization. Notice that A(x, u) has a 
double pole at u = 0, which is converted into a single pole by the factor 
sin7TM coming from the transform of g(a) (cf. eq. (|39|)). The softening of the 
singularity occurs because infrared singularities come from the integration 
over x around the elastic region. The complete box diagram involves instead 
the Borel transform of 

(46) 



given by 



and 



H 2 



u 

\ for {p-kf < (47) 

u 

f or ( p - k) 2 > (48) 

\{p — ky ) 

In this case the factor sin7ru is absent and there is a double pole at u = 
in the final result, coming from the product of the infrared and the collinear 
singularity. 



(p — ky 



K 2 



4 The kernel 

Collecting formulas (p5D , ([34]) and (f44|), we have for the total rate: 

w = w + w a + w b = (49) 
2E ^ q 2 Tr[xp7^(xp + g)7 M ] q(x,Q 2 / /x 2 ;t) 

where q(x, Q 2 / fi 2 ; t) is the parton density of a quark in the Borel plane, 

= w-'WO- ^mff ct£p A(x ' u) - (50) 

The evolution kernel is computed by taking the derivative with respect to 
log Q 2 of the parton density: 

P qq {x,Q 2 /yL 2 -t) = Qi^L q ( Xj Q*/^t) = 



12 



Qn 3 b N \Q 2 J (l-rr)"- 1 " " "v~> - 

Notice that the evolution kernel in the coupling constant space and in the 
Borel space are computed in the same way, because taking the derivative 
with respect to log Q 2 commutes with the Borel transform. 

The factor u sin ttu cancels the single as well as the double pole at u — 
in A(x, u). The evolution kernel is therefore finite at u — 0, as it should 
be for an infrared safe quantity. Notice that the parton density itself has a 
simple pole at u — 0, coming from the collinear singularity. The function 
u sin 7iu has simple zeros at the integers u = l,2,3...,ra, .... The singularities 
of the Borel transform on the real axis originate therefore only from the term 
proportional to the LerchPhi function: 

P„(x,t) = ^iv(^) {l X __ x)u+1 u^u LerchPhi(-a;,2,- M ) 

+ (regular function) (52) 

There are ultraviolet renormalons corresponding to simple poles at the inte- 
gers u — 1, 2, 3, n, .... The leading renormalon is at u — 1 and contributes 
to the evolution kernel by a term of the form: 

2 /j 2 x 3 1 
Pgq ^ = 3n 2 boN Q 2 (l-x) 2 + ^ high6r ° rder renormalons ) ( 53 ) 

The leading non-perturbative correction coming from the ultraviolet renor- 
malon at u — 1 is therefore given by: 

"■«<*• »> = " sip (T^F ^ < 54 » 

where we have taken the pole residue as a measure of the ambiguity. 

We can have a qualitative picture of the ambiguity in QCD by changing 
the sign of the /^-function, i.e. considering the case bo < 0: 

6P„(x,a) = - j^-y 2 ^ e-VNK, 2 ) + ( h .o. corrections) = 

- 3^g_ ^ + (h.o. corrections) (55) 

where by higher order corrections we mean terms of the form (Aq CD / 'Q 2 ) n 
with n > 1. 
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5 Discussion 



Eq. (]55|) is our main result. It says that the evolution kernel has corrections 
of order 1/Q 2 . The most interesting result is that the coefficient of the power 
correction is a function of the Bjorken variable x, and scales as 

(56) 



'1 — x 



The behaviour (^) does depend on the dynamics of the model and cannot be 
derived with kinematical methods. Non-perturbative corrections grow with 
x and diverge as 1/(1 — x) 2 for x — > 1, i.e. in the elastic region. This is a 
stronger divergence than that one of the perturbative kernel (computed up 
to two loops |TIJ), given by 1/(1 — x). There exists therefore a critical value 

at which non-perturbative effects dominate over perturbative dynamics and 
the parton model becomes irrelevant. There is a consequent breakdown of the 
factorization in the cross section (^). We believe therefore that our analysis 
has not only an interest of principle but may have also phenomenological 
applications, identifying the domain of perturbative QCD. The values of x c 
are quite close to one for any reasonable value of Q 2 . For example, x c ~ 0.97 
at Q 2 = 3 GeV 2 , while x c ~ 0.999 at Q 2 = 100 GeV 2 . 

On the contrary, non-perturbative corrections vanish as x 3 for x — > 0. 
They are therefore completely negligible with respect to the perturbative 
terms, which do not vanish in the small x region. This conclusion however 
has to be taken with care because we believe that QED at order 1/N is not 
completely appropriate to describe small x physics. 

The qualitative features of (|56|) can be understood with the following 
qualitative considerations. As x — > 1, the invariant mass of the hadronic 
(partonic) final state and the typical virtuality of the gluon propagator k 2 
reduce progressively. The effective coupling constant a s (k 2 ) is therefore eval- 
uated at a scale progressively closer to the pole: k 2 — > A 2 . The result is that 
non-perturbative effects increase without bound. On the contrary, the invari- 
ant mass of the final state and the typical virtuality of the gluon propagator 
diverge as x — > 0. The coupling constant is evaluated at a scale progressively 
larger with respect to the pole, and non-perturbative effects tend to zero. 
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6 Conclusions 

We analyzed non-perturbative effects in deep inelastic scattering of a lepton 
off a hadron. The idea is that the complete perturbative series hiddens non- 
perturbative properties of the process. Generally speaking, these corrections 
turn out to be small, thereby validating the consistency of the perturbative 
approach. I am also able to explain semiquantitatively such phenomena as 
the precox scaling, i.e. the scaling at quite small momentum transfer. 

There is however a strong divergence of the non-perturbative terms near 
the elastic region. This implies a breakdown of the parton model, as it is 
conventionally understood. 

We believe that our analysis can be extended to other non infrared-safe 
processes, such as 7 — 7 or hadron-hadron collisions. 
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A Borel transform integrals 

In this appendix we describe the computation of the integrals 7i, J 2 and J 3 , 
given in eqs. d43|) . We also prove that the Borel variable u is a regulator of 
the soft divergences (i.e. infrared and collinear). 

For concreteness, let us consider 1\\ the computation of I 2 and J 3 is analogous. 
Making use of the identity 




(58) 



the integral can be written as 



h 




(59) 



where 




(60) 
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J(fi 2 ) is a two body phase space integral with a collinear singularity regulated 
by the 'gluon mass' /i ^ 0. The original integral I\ involves an integration 
over /i 2 with a weight function given by ( / u 2 )~ u_1 . By taking u negative and 
large enough we can suppress the contribution of J(fi 2 ) at small /i 2 , where 
the singularity occurs, in such a way to render the integral finite. Since (as 
well known) the singularities are at most powers of logs, it is sufficient to 
take u < 0. 

J(/i 2 ) is computed by changing the loop momentum to k' = p — k and 
going in the reference frame in which p + q = and q is directed along the 
z axis. In this coordinate system po =\ p | = | q |= q s and q ■ k = —p ■ k = 
q s k cos 6. Integrating the massless delta function, the integral becomes: 

J( 2 x f d 3 k 5 + ((q + q s f - 2(g + q s )k - /i 2 ) 
l/i ' J 2k q 2 - 2q k - 2q s k cos 9 ~~ 1 ' 

Going to an adimensional momentum defined by y = 2k/ \q + q s ) and inte- 
grating the massive 5- functions, we have the result: 

= __1 -log (j_Z^!L (62) 

^' 4q s (q + q s ) S 2 - (1 + v)r 1 ' 

where v = q /q s < 1 and r = fi 2 /(q + q s ) 2 . 

The finite terms for /i 2 — > generate simple poles at w = in 7i, while the 
logarithmic term generates a double pole. 

We perform now the last integration over /i 2 and express go an d q s in 
terms of x and Q 2 according to the formulas 

whose inverse are 



Qs = 





The result is: 

7T X U+2 

h — 



2 (Q 2 )» +1 (l - xY' 1 



1 1 



, -F(l,l- U ,2-u;l-a;) 
x(l — X) U X u[u — 1) 

(65) 
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where F= 2Fi(a, b, c; z) is the hypergeometric function, 
F(a, b, c; z) 



a b a(a + 1) b(b + 1) 9 
1 H tz + 7——— — ——z + ... 



c-1 c(c+l)-l-2 



(66) 



The evaluation of I 2 is analogous to that one of Ii. To simplify the angular 
integration, it is convenient to write 2k ■ p = k 2 — /i 2 . The integral splits as 



where 



and 



h + Io 



, 1; «5 + ((fc + 9) 2 )* + ((p-*) 2 -^ 2 ) 



(67) 
(68) 
(69) 



K(fi 2 ) has a power singularity of the form l//i 2 for /i 2 — > which, integrated 
together with the factor (fi 2 )~ u , generates a pole at u — 0: 



7T 



X 



u+1 



2 (Q 2 )« +1 (l -x)"- 1 
The total result is: 

r U+l 



xu u — 1 



F(l,l-u,2-«;l-a;) 



(70) 



7T 



/ 2 = - 



X" 



2 (Q 2 ) u+1 (l - xY' 1 



'1 — x) V? 



+ 



1 



u + 1 



X u u(u — 1 , 



F(l,l-u,2-«;l-a;) 



(71) 



I2 has the same degree of singularity as I\ at u — 0, given by simple and 
double poles. This implies that k 2 and 2/c • p have to be considered of the 
same order in evaluating the collinear singularities. 
For I 3 we have the formula: 



7T 



X 



u+2 



2 (g 2 )«+ 2 (i - xy- 1 
4 



■2 , l ° gX -u,l-u- -x) 



, 1 — x) 2 u 



[i-xy 



LerchPhi(— x, 2, —u) — 



(i-xy 



;h(x, u) 



(72) 
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where LerchPhi(z, s, a) is the Lerch trascendent of $, defined by: 

oo 

LerchPhi(z, s, a) = > -. — , (73) 

with the terms for which a + n = omitted, and h(x, u) is a function defined 
by the integral: 

Hx,u) = fdy ^-?-*™ . (74) 
v ; 7o [1 + xy] K J 

This function has simple poles at the integers: 

w \ /-, \ 1 (l — x)(l — 3x) 1 

h(x,u) = (1-x) J -\ - + ... 75 

u — 1 2 -u — 2 

B Planar gauge 

In this appendix we prove that collinear singularities decouple from the final 
quark leg in the planar gauge (|24]) , so that the only diagrams contributing 



to P qq are box diagrams |]TT| . In the collinear limit, 

k = xp, (76) 

the polarization sum becomes 

Pfi(xp + q) v + p v (xp + g) M 



Spvip -k,n) = -g^ + 



p-q 

,2 



+ corrections of order k . (77) 



It is a projector orthogonal to the scattering hyperplane, i.e. the hyperplane 
spanned by p and xp + q: 

Pf ,S^(p-k,n) = 0, {xp + q)^S^{p-k,n) = 0. (78) 

The trace coming from the fermion line is a tensor T^ u lying in the scat- 
tering hyperplane in the limit (|76|): T^ v = T^ u {p 1 xp + q). Therefore, the 
contraction with S^ v vanishes: 

T^S^ = at k 2 = 0. (79) 
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The diagrams with a photon radiated by the final quark contain a single 
quasi-real propagator, and the integrals are of the form: 



/rib* 
^ (80) 

Because of eq. (0), these diagrams do not contain collinear logs. 

In the box-diagrams instead, there are two quasi-real quark propagators, and 

the integrals are of the form 

The terms in T'^S^ proportional to k 2 give rise to the collinear log. C.V.D. 
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FIGURE CAPTIONS 

Fig.l: bubble summation in the photon propagator; 

Fig.2: lowest order diagram for the scattering of a quark in an external e.m. 
field; 

Fig.3: diagrams of order a for the scattering of a quark in an external e.m. 
field; 

Fig.4: Diagrams of order 1/N for the scattering of a quark in an external e.m. 
field. 
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